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We use the post-Newtonian- Affine model to assess the validity of the adiabatic approx- 
imation in modeling tidal effects in the phase evolution of compact binary systems. We 
compute the dynamical evolution of the tidal tensor, which we estimate at the 2 PN 
order, and of the quadrupole tensor, finding that their ratio, i.e. the tidal deformabil- 
O^l 1 ity, increases in the last phases of the inspiral. We derive the gravitational wave phase 

corrections due to this phenomenon and quantify how they affect gravitational wave 
detect ability. 
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1. Introduction 



Coalescing binary systems composed of neutron stars (NS) and/or black holes (BH) 
, are one of the most promising sources of gravitational waves (GWs) for second 

and third generation of ground based interferometers [1]. Additionally, certain key 
. features of NS tidal deformations may provide precious information about the NS 

I/*-) \ equation of state (EOS) [2-5]. In current semi-analytical approaches, finite size 

effects are treated by assuming that the NS deformability is described in terms of a 
set of quantities, the Love numbers [6,7], among which the apsidal constant k2 is the 
most important one for GW phenomenology. Within this framework and under the 
so called adiabatic approximation, the evolution of a star embedded in an external 
quadrupolar tidal field, Cy, is governed by the equation 

Qij = -(2/3)fcai& s Ck , (1) 

where Qij and i?NS are the star quadrupole moment and radius at isolation, respec- 
tively. 

In this work, on the basis of the results of Maselli et al. 2012 [8], we employ 
the post-Newtonian- AfBne model (PNA) [9], recently developed to study tidal in- 
teractions in BH-NS and NS-NS systems, to asses the range of validity of Eq. ([1]). 
By computing the ratio ~ Qij/Cij as function of the orbital separation for several 
binary configurations, we find that k^ is not constant during the inspiral and that 
it grows in the last stages of the inspiral. Finally, we estimate the impact of such 
effect on future GWs detections. 

2. The PNA model 



To describe stellar deformations, we improve the PNA model, based on the Affine 
approach [10-14], which relies on the assumption that a NS in a binary system 
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is deformed into an ellipsoid by the tidal forces exerted by its companion. Under 
the Affine hypothesis, the infinite degrees of freedom of the stellar fluid reduce to 
five dynamical variables, namely the three principal axes of the ellipsoid and two 
angles related to the star angular velocity and the internal fluid motion. The main 
improvements to the original Affine approach [10-12], were brought in [9,13,14], 
where the NS spherical equilibrium configuration was determined by solving the 
relativistic stellar structure equations, and a more careful treatment of coordinate 
systems and of relativistic corrections to the star internal dynamics were introduced. 
Furthermore, starting from a two-body post-Newtonian metric [15,16], an explicit 
expression for the tidal field was derived. This represents a major change from the 
original approach in which the NS moves in a background space-time described by 
the Kerr metric, and it allows to consistently treat stellar deformations in NS-NS 
binaries. The orbital motion is handled within a post-Newtonian framework, as 
well: the phase evolution is determined by means of the TaylorT4 approximant [17], 
including the effects of tidal interactions on the orbital motion, up to 1 PN order [18]. 

The source of stellar deformations in the PNA model is given by the quadrupolar 
tidal tensor CV, = Ra^se^e^eJ^e^y where R a {s-ys is the Riemann tensor of the 
3 PN metric describing the orbital motion, and e(j = o....,3) is an orthonormal tetrad 
field fixed to the star center of mass and parallel transported along its motion. The 
results obtained in [9], where wc derived CV, up to C(l/c 3 ), were extended in [8] 
to include C(l/c 4 ) terms. With this aim, we also computed the 2PN contribution 
to the spatial tetrad vectors e| fc N introduced in [19]. Cy depends on the masses of 
both compact objects mi, 2, and includes, 1.5 PN order, the spin contribution from 
rotating objects. 

3. Tidal Love number evolution 

We solve the PNA equations of motion up to the onset of the NS mass shedding, for 
a representative set of binary system configurations, by considering two equations 
of state which are expected to cover the range of possible NS deformabilitics. For 
each binary model, we compute ki as the ratio between the quadrupole and the tidal 
tensor, from Eq. fl}. We find that this quantity is a function of the orbital separation 
r, i.e. &2 = ^2(^)5 that increases during the binary inspiral by a ~ 10% — 30% 
factor. The constant Love number ki should be regarded as the asymptotic limit 
of this function, i.e. &2 = lini r ->-co &2 (r) ■ This effect is strongly dependent on the 
EOS choice: it is larger/smaller for a less/more compact NS. The dependence of 
this effect on the mass ratio is weaker, and it increases more for larger values of 
q = 7712/ mi < 1- This effect appears to be nearly insensitive to the NS mass. 

Replacing the Love function fc 2 (r) into the expression of the GW signal, we 
estimate a new correction to the GW phase. This adds linearly to the point-particle 
(PP) term, and to the tidal term (T) computed by assuming the Love number to 
be constant: 



(2) 
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where 5ipT describes the growth of k%{r) during the inspiral. We find that the new 
contribution is always less than 1 radian for all considered configurations, up to the 
onset of mass shedding. 

To assess the validity of the adiabatic approximation as far GW detectability is 
concerned, we compute overlaps [20] between PP templates and real signals which 
take into account tidal effects at the best of our knowledge, i.e. including Sipr- We 
find that in the worst case, i.e. for a NS-NS binary with a standard total mass of 
2.8M©, the number of missed events due to the use of PP templates could be as 
high as 11% and 13% for the second and third generation of ground based detectors, 
respectively, depending on the NS equation of state. By using templates which in- 
clude tidal corrections by means of the constant Love number, however, the overlaps 
would differ from unity by less than a part in one thousand for all configurations 
considered. We may thus conclude that building gravitational waveforms within the 
adiabatic approximation is very reliable up to the mass-shedding. 

4. Conclusions 

In this work we reviewed the results of [8], in which we improved the PNA model 
by computing the 2 PN contribution to the post-Newtonian tidal tensor. We used 
this framework to assess the reliability of the Love number adiabatic approximation 
and found that the NS deformability increases during the last phases of the inspiral 
up to 30%, depending on the NS EOS. 
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